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^ . Abstract 
^ ■ 

Q_i' In this paper we present a calculation of the AS = 1 effective weak 

Hamiltonian including next-to-leading order QCD and QED correc- 
D i tions. At a scale /U of the order of few GeV, the Wilson coefficients of 

' the operators are given in terms of the renormalization group evolu- 

■ tion matrix and of the coefficients computed at a large scale ~ My/. 
^ , The expression of the evolution matrix is derived from the two-loop 

^ \ anomalous dimension matrix which governs the mixing of the relevant 

current-current and penguin operators, renormalized in some given 
regularization scheme. We have computed the anomalous dimension 
matrix up to and including order and aeCtg in two different renor- 
malization schemes, NDR and HV, with consistent results. We give 
many details on the calculation of the anomalous dimension matrix 
at two loops, on the determination of the Wilson coefficients at the 
scale Myi/ and of their evolution from Mw to ^. We also discuss the 
dependence of the Wilson coefficients/operators on the regularization 
scheme. 



> 



^Unite Propre du Centre National de la Recherche Scientifique, associee a I'Ecole Nor- 
male Superieure et a rUniversite de Paris-Sud. 



1 



1 Introduction 



In this paper we present a calculation of the two-loop anomalous dimension 
matrix relevant for AS = 1 decays. The anomalous dimension matrix in- 
cludes leading and sub- leading corrections at orders a^, a^, and a^ae- 
The calculation has been performed in two different regularization schemes: 
naive dimensional regularization (NDR) and 't Hooft-Veltman regularization 
(HV)!!!]. We verify that the results obtained in the two schemes are com- 
patible both in the strong and electro-magnetic case. We give many details 
on the calculation itself, on the definition of the renormalized operators, on 
the relation between different regularization schemes, on the role of counter- 
terms for operators which vanish by the equations of motion, on the gauge 
invariance of the final result etc. We also report a list containing the double 
and single pole contribution of all the diagrams, counter-terms and "efferves- 
cent" operator counter-terms in both the schemes used in this calculation. 
The list of all the diagrams may be useful to check our results and for further 
applications. 

A comparison of our results with a parallel calculation reported in refs. 0- 
[^] is also presented. We agree with the authors of these references in the 
NDR scheme. For some diagrams computed in the HV regularization scheme 
we however disagree and we explain the origin of the difference. In the electro- 
magnetic case we show that, by using the values of ref. for the diagrams 
computed in HV, it is not possible to satisfy the expected relation between 
the two-loop anomalous dimension matrix computed in NDR and HV. On 
the contrary our results satisfy the expected relations for both terms of order 

and asCte- We have also checked the consistency of our calculation at 
order a^. 

The authors of ref. have only computed the anomalous dimension in 
NDR and derived the result in HV by using the one-loop anomalous dimen- 
sion and coefficient matrix. We thus agree with their final result for the 
two-loop anomalous dimension matrix in HV too, in spite of the different 
results for some diagrams between this work and ref. 0. 

In this paper we have preferred to report only the calculation of the 
two-loop anomalous dimension matrix with as many details as possible and 
postpone a discussion of the numerical calculation and uncertainties for the 
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coefficient functionsQto a separate publication 0. A phenomenological analy- 
sis of e'/ e, using the results reported in this paper, has already been presented 
in ref.0 and will be extended/upgraded in ref.]^. 

The paper is organized as follows. In sec.Q we introduce the general 
formalism. The calculation of the coefficient functions at the W scale is 
summarized in sec.^. A detailed discussion of the scheme dependence of the 
anomalous dimension matrix and of the relation between different renormal- 
ization schemes is given in sec.^ We also present a convenient definition of 
the renormalized operators which makes the evolution matrix scheme inde- 
pendent. This definition may be useful to predict weak amplitudes by com- 
bining the Wilson coefficients with the matrix elements of the corresponding 
operators computed with a non-perturbative method, as for example lattice 
QCD|^, 1^. In the same section we also recall some basic features of the HV 
and NDR regularizations. Sec.^ is the main section of this paper. There 
we describe the calculation of one- and two-loop diagrams, discuss the role 
of the so called "effervescent" operators and comment on the subtraction of 
counter-terms corresponding to operators which vanish by the equations of 
motion. Double and single pole contributions of all the relevant Feynman di- 
agrams, in the NDR and HV regularizations are given in the Appendix. From 
the calculation of the Feynman diagrams, the two-loop anomalous dimension 
matrix is derived and given in the NDR and HV schemes in sec.|^. 



2 General Formalism 

Effective Hamiltonians for non-leptonic decays of hadrons composed by light 
quarks {K, D and B mesons for example) are defined by Wilson operator 
expansions of products of weak currents|0]-0: 

< F\neff\I > = g'w/8 J d^xDw{x\ M^) < F\T Jt(0)) |/ > 

^ 5:<F|g,(^)|/>Q(/i) (1) 

i 

For kaon decays, in the limit in which we neglect quark masses, only four- 
quark operators appear on the r.h.s. of eq.(|I|). The AS* = 1 effective hamil- 

^ Dependence of the coefBcient functions on the renormahzation scale /i, on Aqcd and 
on the renormahzation prescription for example. 
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tonian can then be written as : 

Kfr' = A„^[(l-r)(Ci(/i)(Qi(/x)-Q^(/i)) + C2(/i)(Q2(/^)-Q2(/^)) 
+ TY,Q,{^^)C,{^^)\ (2) 

i=l 

where A„ = VudV*^ and similarly we can define Ac and A^. r = —Xt/Xu and 
Vij is one of the elements of the CKM|TO|, pJ]] mixing matrix. The operator 
basis is given by: 



Qi = 


(SqM/j) (y-A) {Uiida)(y^A) 


Q2 = 


{SaUa)(V-A) {u 13d ji) (y - A) 


Q3,5 = 


{Sada)(V-A) ^ {Wp) (V^fA) 
q=u,d,s,--- 


Q4,6 = 


{Sadp)(v-A) X! i.Wa){VTA) 
q=u,d,s,--- 


Q7,9 = 


3 

-{Sada){V~A) X! ^qi^lilp) {V±A) 
q=u,d,s,--- 


Qs.io = 


3 

■^{Sadf3){V-A) ^q{Q!iQa){V±A) 
q=u,d,s,--- 


QI = 


(Sq-C/j) (y „A) {Cl3da){V-A) 


= 


iSaCa){V-A) (C/jrf/j) (V-A) 



(3) 

when QCD and QED corrections are taken into account p|, 0, [|l2|-[llll- In (H) 
the subscript {V ± A) indicates the chiral structure and a and j3 are colour 
indices. The sum is intended over those flavours which are active at the scale 
fi. We have completely ignored the effects due to the operators Q\ and Q2 
which are the analog of Ql and Q2 with the charm quark replaced by the 
bottom quark. In ref. ||12| it was indeed shown that these operators have a 
negligible effect on the evolution of the Wilson coefficients of the operators 
@. Their inclusion, once that the anomalous dimension matrix is known, is 
in any case elementary. For rric < ^ < nib we have used the relation: 

Qio = Qg + Q4 - (4) 

to eliminate Qio from the evolution equations. We think that expression (|^) 
is the most trasparent for fi > rric. It shows that we can find all the Wilson 
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coefficients by evolving Ci^io (C*i,9) via a 10 x 10 (9 x 9) evolution matrix down 
to fi = rUb {jric < /i < rrih)- The generalization to fi < rric is straightforward 
and can be found, for example, in a recent paper on e'/e|TB[. 



The operators Qi{n) are renormalized at the scale n in some given scheme. 
The corresponding coefficients Cj(/i) are scheme dependent. The dependence 
on the regularization scheme appears at one loop, when we express the orig- 
inal current-current product in terms of the Wilson operator product expan- 
sion (OPE), see eq.®. 

To obtain the coefficients Cj(/i) at next-to-leading order (NLO) two steps 
are necessary: 

1) The calculation of the coefficients at a given scale, for example M\y 
or rrit, including corrections of order and Oe- 

2) The calculation of the two-loop anomalous dimension up to 0{a^) and 
O(a,ae)0- 

The results of steps 1) and 2) depend on the regularization scheme and 
on the normalization conditions imposed on the renormalized operators, as 
will be discussed below. We have done our calculations in two popular reg- 
ularization schemes, i.e. the t'Hooft-Veltman (HV)|]ll and the naive (NDR) 
dimensional regularization schemes. In both cases we have obtained the 
renormalized operators via the standard modified minimal subtraction pro- 
cedure MS. We will also discuss other renormalization prescriptions which 
can make the renormalization group evolution matrix scheme independent. 

In presence of 75 and in certain regularization schemes, the axial vec- 
tor current may develop an anomalous dimension at the two-loop level. In 
defining the evolution matrix one has to take into account this effect. Al- 
ternatively one can impose to the current a certain one-loop renormalization 
condition such that its two-loop anomalous dimension is zero. We prefer this 
second solution and discuss this point in sec.0. 

To make easier a comparison with previous calculations on the same 
subject, 0-0] and [0, we follow as close as possible the notation introduced 
in ref.0 and write: 

n^ff'-Q^{f^)cif^) (5) 



Through this paper we neglect terms of order 
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where is a row vector whose components are the operators Qi.io of the 

basis and C{fi) is a column vector, whose components are the correspond- 
ing Wilson coefficients. C{ii) are expressed in terms of C{M]y) through the 
renormalization group evolution matrix iy[/i,Mi4/]: 

C{fi) = W[fi,Mw]C{Mw) (6) 
The coefficients C(/i) obey the renormalization group equations: 

(-| + f^i^s)^ - ^^^^^)Cit, aM, «e) = (7) 

where t = ln{Myy/fi'^) and we ignore the running of ae- The factor of 2 in 
eq.(|^) normalizes the anomalous dimension matrix as in refs-^-jQ. 

7 = 7q - 2 7j i (8) 

is the anomalous dimension matrix of the operators minus twice the anoma- 
lous dimension of the weak current in a given renormalization scheme. In 
eq.(H) 1 is the identity matrix. 

To simplify the discussion, we first consider the case where there is no 
crossing of a quark threshold when going from Myy to fi. The relevant for- 
mulae for the general case will be given at the end of this section. At the 
next-to-leading order, by expanding Mvi/], we can write: 

W[H, Mw] = M[fi]U[ij, Mw]M'[Mw] (9) 

with: 

MM^(l.^/<-)(i.^i)(i.^p) (10, 

and 

M'[M^] = (l- -^P] ( i - ^^iiiM A (i _ ^k) (11) 



as{Mw) J \ 47r y V 47r 

We substitute the expression of C(/i) given in eq.(D in the renormal- 
ization group equations (0), using iy[/i,Mt^] written as in eqs.(P-|TT|). By 
expanding the anomalous dimension matrix, which includes gluon and pho- 
ton corrections, up to order and aea^: 

7 = ?^7f + + i^rii'^ + ^^li'^ (12) 

47r Att Mtt^ A-kA-k " ^ ^ 
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we obtain the expression for U : 



U[lJ,,Mw] = Ta,exp(- f 

J a 



(o)r 



as{Mw) 



(13) 



in the basis where 7^°^ is diagonal. T^^ is the ordered product, with increasing 
couphngs from right to left. The matrices P, J and K are solutions of the 
equations^: 

(14) 



P 



P 



7, 



Wo 



J 



J. 



1. 



(0)T 



2/?o 



/3: 



(1)T 



'2 7s 



2/?o 



(15) 



{o)r 
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«^ + 7f^^ + 7i^)^P+f7f^,iP 



2/3iP-^P7f^ (16) 



In eqs. ([T^-[T6|), /5o and /5i are the first two coefficients of the /3-function 
of as- U and P are determined by the leading logarithmic (LO) anoma- 
lous dimension matrices 7^°-* and ^f"* and are regularization scheme indepen- 
dent. On the other hand the two-loop anomalous dimensions 7^^-' and 7^^-', 
and consequently J, K and W[ii, Mw], are regularization scheme dependent. 
Eqs.( p!lHT^ ) can be easily solved in the basis where 7^°^ is diagonal. The so- 
lutions develop singularities which however cancel in the final expression of 
W[ii, Mw\- It is indeed possible to find an explicit form of Mw] which 
is not singular. This form was used in the numerical calculation of ref.|@]. 

The initial conditions for the evolution equations, C{Mw) are obtained 
by matching the full theory, which includes propagating W , and six 
quarks, to the effective theory, where the W , Z° and top quark have been 
removed simultaneously. In general, C{Mw) depend on the definition of the 
renormalized operators in a given regularization scheme. A scheme indepen- 
dent way of defining C{Mw) and the evolution matrix will be discussed in 

^ The last term in eq.(15) of ref.|^, corresponding to eq.(p^ of this paper, contains a 
misprint. The expression which is reported here has however been used in all the numerical 
calculations of the Wilson coefficients of ref.|^, which are consequentely correct. We thank 
A. Bur as for finding the misprint. 
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sec.^. In the full theory, all the coefficients coming from current-current and 
penguin operators have been computed in refs.[12|,|20|-|22|. In this work we 



have only computed, in the effective theory for the HV and NDR regulariza- 
tion schemes, the 0{as) and 0{ae) one- loop corrections necessary to impose 
the matching conditions on C(Miy). 

We now discuss the modifications to the evolution matrix in presence of 
a heavy quark threshold. These modifications are necessary when /i < mb- 
The estension to the case when n < rric is straightforward. The matrix 
W^[/ii,/i2] depends on the number of active flavours / in the interval [fii,fi2]- 
We denote it by Wf[fii,fi2]- On the other hand when we cross a threshold, 
the renormalization conditions of the operators are in general changed^. At 
threshold it is thus necessary to introduce a suitable matrix T 0. This matrix 
allows for the matching of the evolution between scales larger and smaller 
than the threshold. Thus for example, to evolve the coefficients from Mw to 
< nib, one has to use the following expression [p^j : 



W[^i, Mw] = W4fi, mfe] WsK, Mw] (17) 

where 

f = i + ?^Sf^ + ^6s^ (18) 



The matrices 5f and 6s relevant for eq.(PTD are given in sec.^TT] together with 
all other one-loop results. 



3 Calculation of the Coefficients C{M] 



w, 



In order to compute the Wilson coefficients of the OPE at a scale yU ~ M\y 
{rrit), we have to consider the full set of current-current, box and strong, 
electro-magnetic and penguin diagrams up to and including O^Gpc^s) 
and 0(G'i?ae)[|Ill!llO|-|2^. In the current-current case we compute the di- 



agrams with external momenta \pf\ ^ fi^ (i.e. we neglect terms of 

order f^l/M^r) with massless external quark states. The dependence on the 

^This however does not happen with scheme independent renormahzation conditions. 
^ The matrix T is different for different thresholds, depending on the quarlc electric 
charge. 
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external momenta only appears in logarithms, proportional to the anoma- 
lous dimension of the operators, and in the operator matrix elements. Strong 
and electro-magnetic penguin diagrams are also computed with external mo- 
menta \pj\ <^ (Ip^I ^ m^) and massless external states. In the case of 
penguin diagrams, the logarithmic dependence on the external states appears 
as a dependence on the momentum transferred through the gluon or photon 
propagators, ^ fJ^l- 2'°-penguin and box diagrams can be computed with 
zero external momenta and including only the top quark contribution, since 
they are infrared finite, cf. B{xt) and C{xt) in table |l|. 

We now introduce the notation necessary for the calculation of the co- 
efficients C{Mw). In the full theory, the direct calculation of the current- 
current, box and penguin diagrams at one loop (including order and 
corrections) has the form: 

< JJ > ~ < Q(°)^ > ff + + ^D(^)" 

L 4v]- 47r J 

= < Q^{Mw) > C{Mw) (19) 

where < Q^^^^ > are the tree-level matrix elements and T^^^ and D^^^ de- 
pend on the regularization scheme and on the external quark states. Indeed 
all the W-g and W-'j box diagrams, being finite, are regularization scheme 
independent. The axial vector vertex diagram however does depend on the 
regularization scheme. The diagrams necessary to obtain T^o)'^^) and D^^^ 
are shown in figs. 1,2 and 4-8. By inserting the renormahzed operators of the 
effective hamiltonian (|]) in the diagrams reported in figs. 3 and 9, we then 
compute the one-loop current-current, strong and electromagnetic penguin 
diagrams between the same external states, using the same regularization 
scheme. In this case we obtain: 



<Q(Mh/)>= l + -^f + -^S <QW> (20) 
^ 47r 47r ' 

T^^^ and D'^^^ contain logarithms of the external momenta, whose coeffi- 
cients are proportional to the one-loop anomalous dimensions of the relevant 
operators (the W and top masses acting as an effective ultraviolet cutoff). 
Thus for example T^^^ contains terms ~ ln{My//n'^) (~ ln{m^ / fi^)) . On 
the other hand the insertion of the renormahzed operators Q{^) between the 
same external states goes like ln{fi'^/ /Xq), where /i is the renormalization scale. 
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If we choose the renormahzation scale /i = Mw, the logarithms disappear 
when we compute the coefficients C{Mw), which are obtained by comparing 
eq.([T|) with eq.(|2g): 

C{Mw) = f (0) + ^(f - f^f + ^fD(i) - s^f (21) 

ATT ^ ' 4:71 ^ ' 

We take as a fixed coupling constant and in eq. ( |2TD has to be interpreted 
as as{Mw)- T^^'^ and (D^^^ and s^) depend on the external states and 
on the regularization scheme. However their difference depends only on the 
regularization scheme. We will give the results for both the HV and NDR 
regularizations. 

We now describe separately the strong and electro-magnetic current- 
current, penguin and box diagrams. 



3.1 Current- Current 0{as) Diagrams 

The coefficients C{Mw) depend on the combination T^^^ — f^T^^^ which is 
independent of the external states. Thus we can choose to compute the 
diagrams in figs. 2a-2c and 3a-3c with external states different from diagram 
to diagram (but equal for corresponding diagrams in the full and effective 
theory, 2a and 3a for example). We have choosen the external momenta as 
shown in the figures [|. At one loop in a^, < J J > can only mix with the 
operators Q12 and Ql 2 of the list given above. We only discuss the mixing 
of Qi,2, since the case Ql 2 proceeds in the same way. 

When we compute the vertex corrections to the weak charged current, 
fig. 2a, and combine it with the renormahzation of the external quark lines, 
fig. 4, the axial current is in general subject to a finite renormahzation which 
depends on the regularization scheme. A finite correction to the axial current 
at one loop implies that the current have a non-zero two-loop anomalous 
dimension. The current anomalous dimension must be subtracted from the 
anomalous dimension of the operators of the effective Hamiltonian, see eq.(H). 
Alternatively we can apply a finite renormahzation to the current in such a 
way that its two-loop anomalous dimension is zero. This procedure modifies 

^ The same choice of external momenta for t he c urrent-current diagrams has been 



adopted for the 0{ae) corrections described in secj3.2 
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the coefficients C{Mw)- It is equivalent to impose that the current obeys the 
Ward identity which states that, in perturbation theory, the axial current 
must be conserved in the limit in which the external quark masses vanish. 
The two definitions of the axial current will give, at the NLO, the same 
physical result. We prefer to choose the second alternative, i.e. 7j = at 
two -loops, and avoid the subtraction of the current anomalous dimension. 

With our choice of the renormalized weak current, in the HV scheme, we 
find: 

'^HV-ccgWlw) - — 

^(2) (M \- ^ (00) 



where the subscript ccg indicates that this is the contribution from current- 
current diagrams, "cc", from the exchange of a gluon, "^f". In the calculation 
of the coefficient functions (7(Mvi/), as{Mw) has to be understood as the 
running coupling constant, computed in the 5-fiavour theory at the scale M^. 
This is the correct procedure when rrit is larger than My/, as suggested by 
the mass lower bound obtained for a top quark which decays in the standard 
modes. In the NDR scheme, we find instead: 

r« (M ^ - "^sjMw) 11 

'^NDR-ccgV^'^w) - — 



^(2) (. r ^ _ as{Mw) 11 , . 



3.2 Current- Current 0{ae) Diagrams 

The electro-magnetic corrections are more complicated due to the presence 
of the non-abelian diagrams of figs.Sd and 6b and the W-Z^ box diagrams. 
The sum of the vertex and self-energy diagrams, figs. 5a and 6, contrary to 
the 0{as) case, is neither finite nor gauge invariant. This is related to the 
fact that the sum of all the diagrams in figs. 5 and 6 contributes both to the 
renormalization of the Fermi constant Gp, which is universal for quarks and 
leptons, and to the electromagnetic corrections to the effective Hamiltonian 
at order aep4|, This means that we can reabsorb a part of the 0{ae) 
corrections by a suitable redefinition of Gp- We proceed following ref . [p5|] . 
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In the vertex and self-energy diagrams in figs. 5a and 6a, we write the photon 
contribution using the identity: 



^M^(M^-fc2)-i + (fc2-M^)-i (24) 



The contributions coming from the second term in eg. p^) are divergent. The 
divergent part of these diagrams, combined with the divergent terms from 
diagrams in figs.Sd and 6b, gives a universal contribution (equal for quarks 
and leptons) to the renormalization of the weak coupling constant gw and 
W mass. Thus this term can be reabsorbed in the definition of the physical 
Gp SlS measured in /i-decaysp5[. The first term in eq.(0) is cut-off by the 
additional converg ent factor M^{M^-k^)-^. For this reason it can only give 
a finite contribution ~ Zn(M^//iQ), similar to the contribution of the 'j-W 
box diagrams of figs.5b-c. The sum of the terms ~ /n(M^//iQ) gives indeed 
the current-current contribution to 7^°-*. The remaining genuine 0{GFO;e) 
corrections (i.e. not containing ln{M^/ nl) ) are all proportional to Q2 both 
in the effective and in the full theory [^] -p5||. They give thus a correction 
of order ae/47r ~ 10~^ to a coefficient of 0(1) and can be safely ignored. In 
ref.|]13[ they have instead computed the finite 0{ae) corrections by analogy 



with the 0{as) case, i.e. by considering the differences between the diagrams 
of fig.3a-c and 5a-c with a photon exchanged. In this case one would have 
found: 

C'^'un-cc^iM^) = (25) 

We observe that the above coefficients do not exaust the 0{ae) corrections. 
For example they do not contain the terms coming from the Z^-W box dia- 
grams which exist only in the full theory. 



3.3 QCD Penguin Diagrams 



Penguin diagrams have been extensively studied in the literature |^-|]12l. In 



the full theory, for fiQ ^ nic, using the unitarity of the CKM mixing matrix, 
we have only to consider the contribution coming from the difference between 
the top-penguin and the up-penguin diagrams. The analogous contribution 
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Inami-Lim Functions 


Expression 


B{xt) 


i [xt/{l - Xt) + Xt \\lXt/{Xt - If] 


C{x,) 


\xt [{xt - Q)/{xt - 1) + (3xt + 2)/{xt - If \xixt] 


D{xt) 


-flnxt + (-19x? + 25x?)/[36(xt - 1)^] + 
[xlibxl - 2xt - 6)]/[18(xt - 1)^] \nxt 


E{xt) 


-|lnxt + [x2(15- 16xt + 4x2)]/[6(l-xt)4] \nxt+ 
\xt{l^-llxt-x1)]/[l2{l-xtf] 



Table 1: Basic functions governing the m^-dependence of various weak am- 
plitudes 



from the charm-up penguins is cancelled by the GIM mechanism. We can 
consider the up-quark as massless in the calculation. 

From the diagrams in fig. 8a and 9, in the HV scheme, we find: 

Cfy.,,{M^) = Cfy_^^{M^) = (26) 

where Xt = mj/M^. The modified Inami-Lim function E{xt) is given in 
table ^ In the NDR scheme one finds: 

Cj^ln-.^iMw) = C^nIr-^^A^w) = -^4^ (^^"^^^ ~ D ^^^^ 

These results were already presented in ref.[0]. 
3.4 QED Penguin Diagrams 

In the case of electro-magnetic penguins we have also to consider, besides the 
diagram of fig. 8a, the non-abelian diagram given in fig. 8b. In HV, we find 



13 



the following global contribution for QED-penguin diagrams: 

cSl-miMw) = ^D{x,) (28) 



where the Inami-Lim function D{xt) is given in table |^. In NDR one finds: 



C^NlR-,,iMw) = ^{d{x,)-^] (29) 
6n\ ^ ' 9. 



3.5 Z*^ Penguin Diagrams 

penguin diagrams give 0{ae) contributions] 20 1 as it was the case for 



photon-penguins. The corresponding modified Inami-Lim function van- 
ishes a.s Xt 0. We can then reasonably neglect the up and charm quark 
contributions, even though the GIM mechanism is not active in this case. 
One gets: 



CPv^MM^) = (30) 
^ OTT sm Vw 



cfv_,z«{Mw) = ^iC(x,) (31) 



and the same in NDR. C{xt) is given in table 
3.6 Box Diagrams 

In the approximation of massless light quarks, the box diagrams of fig. 7 give: 



Cg(MH.) = ^-i—2Bix,) (32) 
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both in HV and in NDR. 



A summary of the rcsuhs discussed in this section for the HV and NDR 
regularizations is reported below. We have also included the tree-level con- 
tribution from the diagram of fig.l and neglected terms of 0(Q;e/47r) in C2. 



i) HV 
C2{Mw 

Ci{Mw 
Cio{Mw 

and: 

ii) NDR 
Ci{Mw) = 

C2{Mw) = 



af){Mw) 7 
An 2 

47r 6 



247r 



67r sin -dw 



[2B{xt) + C{xt)] 



af\Mw) 



E{xt) 



E{xt) 



24:71 

Stt 

^ [4C(x,) + D{x,)\ 


Gtt 




(33) 



4C{xt) + D{xt) + -^{lQB{xt) - 4C{xt)) 



sin '&W 



af\Mw) II 



1 - 



47r 2 
af\Mw) II 
An 6 



247r 



OTT 



{em -I 



1 



67r sin 'dw 



[2B{xt) + C{xt)] 
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Cs{Mw) 

CgiMw) 

CioiMw) 



(M, 



247r 

w] 



Eixt) 



E{xt) 





67r 




Stt 

4C(a;t) + {D{xt) 



(34) 



9 
4 

9^ ' sin^d 



1 



105(xO - AC{xt) 



4 Anomalous Dimensions at One and Two 
Loops 

In this section we introduce the notation necessary for the calculation of the 
anomalous dimension matrix 7(05, Oe) in dimensional regularization, eq.(P), 
and recall the rules for the HV and NDR schemes. 



4.1 General Definitions and Scheme Dependence 

For simplicity, we start by considering only one- and two-loop corrections due 
to strong interactions. The modifications, necessary to include the electro- 
magnetic corrections will be given in sec.^ 

The anomalous dimension matrix for the operators appearing in the ef- 
fective Hamiltonian is defined from the operator renormalization matrix: 

7Q(a.) = 2Z-V'^^ (35) 

where Z is defined by the relation: 

Q = Z-'Q"" (36) 
which gives the renormalized operators in terms of the bare ones. 



16 



7Q = ?^7f + 7^7i^^ + --- (38) 



In a dimensional regularization, as in the HV and NDR schemes, from 
eq. (|35|) , we obtain: 

7Q = 2 Z-^ (-ea, + -^Z (37) 

By writing 7q and Z as series in the strong couphng constant^: 

47r ' (47r)2 

^ = 1 + ^Z(^) + 7^^(') + ■ • • (39) 
47r (47r)2 ^ ^ 

we derive the following relations: 

7f = -2eZ(i) (40) 

and 

7(1) = -4eZ(2) - 2/3o^('' + 2eZ(i'z(i) (41) 

where e = (4 — D)/2. Pq is the one-loop coefficient of the /3-function /3(as) 
which governs the evolution of the effective coupling constant: 

(42) 



and: 

An '■^'(47r)2 



(3{as) = -f3o^ -Pi7^ + 0{at) (43) 



Po and Pi are given by: 

_ (lliV-2/) 
Po — 



. ?l^^.fMf-i^f (44) 



^ Since we normalize the weak currennt in such a way that 7j = 0, see below, we have 
7 = 7Q- 
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where / is the number of flavours. The running couphng constant, solution 
of eq.(|D, is : 

_ 1 / f3MH^^''/^QCD)] ^ 

The above equation defines Aqcd at the NLO. 

We can expand Z*^*) in eqs.(^) and (il) in inverse powers of e: 



Z« = E (-)' (46) 



j=0 



The anomalous dimension is finite as e — > 0. This implies a relation between 
the one- and two- loop coefficients of Z\ 

AZf^ + 2/?o^l^^ - 2Zf ^ = (47) 

From the above equations we finally obtain: 

7^ = -2Z« (48) 

and 

if^ = -4Zf ) - 2P,Z^'^ + 2{Z?z'i^ + z'i^Z?) (49) 

We thus conclude that it is sufficient to compute the pole and finite part of 
and the single pole of Z^'^^ in order to obtain the two-loop anomalous 
dimension. Eq.(^) tells us how to derive 7'^^''. In dimensional regularizations, 
such as HV, NDR or DRED (dimensional reduction) however, the calculation 
is complicated by the presence of the so called "effervescent" operators (EO), 
which appear in the intermediate steps of the calculation [1^, |1^]. The EO 



are independent operators which are present in D-dimensions but disappear 
in the physical basis of operators in 4-dimensions. Because of the presence 
of the EO, the products of the matrices zj*^ in eq. (|49|) have to be done by 
summming indices over the full set of operators, including the EO. Only at 
the end of the calculation we can restrict the set of operators to the operators 
of the physical 4-dimensional basis (^). 

Different renormalization prescriptions will define different renormalized 
operators. This happens for example if we adopt the MS subtraction pro- 
cedure in the HV and NDR schemes. Let us denote as Zhv {Zndr) the 
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renormalization matrices in the two cases and thv {^ndr), cf. eg. pO|) , the 
one-loop matrices of the effective theory in the two regularizations. Then we 
must have: 

ZHV = ZNDR{i + ^Af) (50) 

with Af = fj^DR — fuv- Using eq.(^5|), expanded in powers of ag as in 
eqs.(|8D-(|9]), we then find[|: 



Eq.(|5l|) imphes that the combination: 



(51) 



G = 7. 



(1) _ 



(0) 



(52) 



is regularization independent. It can be shown that a consequence of eq. 



is the scheme independence of the combination R 



+ J{S = s^ + Km 



the case of electromagnetic corrections) on the regularization scheme, see for 
example ref. 0. R and S are precisely the combinations which appear in the 
final expressions of the coefficient functions^. 

Using eqs.(PD-(Pl) and ([2l|) , we obtain: 



Cifi) = M[ij]U[fi, Mw]N'[Mw]C'iMw) 
where M[fi] has been defined in eq.(|TO|) and: 



(53) 



w 



a. (Ml 



-P 



w 



An 



Of 



^ ^ A-n Att 



(54) 
(55) 



In the above equation we have neglected higher order terms in or ctg- 
From eq . (|5^) , we conclude that the matrix iV'[Miy] is independent of the 
regularization. We can obtain scheme independent coefficients C(/i) by a 
suitable redefinition of the renormalized operators Q^{fi): 



An 



a 



e 



An 



(56) 



^ The corresponding formula for ji^"^ is Aji^^ 

10 



Af,7, 



(0) 



As, 7, 



(0) 



This is strictly true if we include the terms of 0{ae) of eq.pa). 
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With the above redefinition, one has: 



^ m - (i + ^(/^ + ^^)) (i + ^(i + ^^)) (i + 

(57) 

SO that C(yu) = A^[yu]f/[/i, Miy]A^'[Miy]C"(Mvi/) is regularization scheme in- 
dependent. 

It is not difficult to understand how it is possible to find renormalization 
conditions which do not depend on the scheme. Let us fix the renormal- 
ization conditions by imposing that the matrix elements of the renormalized 
operators have a given value for a certain set of external quark (gluon) states: 

< (3\V\a >= 1 (58) 

This procedure defines the same renormalized operators, i.e. the same co- 
efficients C{fi), in all the regularization schemes^. Notice however that the 
coefficient functions depend now on the external states chosen to fix the 
renormalization conditions, \a > and |/3 > in eg. (|58|) . The external states 
have thus to be specified if one wants to use renormalization scheme indepen- 
dent coefficients. The values of the coefficients can change in a substantial 
way by going from the MS HV or NDR prescription to the scheme inde- 
pendent one. For example, using MS HV we found that Cq is reduced by 
the inclusion of the NLO corrections [^] while it is enhanced in the scheme 
independent case|0, 0. We remark that a consistent treatment of the Wilson 
coefficients and renormalized operators is necessary in order to get, up to 
higher order corrections, the physical result. Such a treatment is possible, at 
least in principle, in lattice QCD but not in other approaches as for example 
the 1/N expansion. 

Expression (|^), which allows us to compute the two- loop anomalous 
dimension in terms of the one- and two- loop renormalization matrices, is 
indeed valid diagram by diagram. This means that we can define the contri- 
bution of any given two-loop diagram to the anomalous dimension by com- 
bining: 

""^^ If the renormalized strong coupling constant differs in two different renormalization 
schemes, as it is the case for HV and DRED for example, the coefficients will have the 
same expression only when given in terms of the same renormalized as 
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i) the single pole obtained from the diagram where we insert a given bare 
operator; 

ii) the single pole of the diagram obtained by substituting to any di- 
vergent sub-diagram the appropriate counter-terms, including those propor- 
tional to "effervescent" operators; 

iii) the single pole coming from the substitution, in any divergent sub- 
diagram, of an appropriate combination of effervescent operators to take into 
account the term 2{z[^^ Z^^^ + Z^^^z[^^) in eq.(||). 

We will call the contribution to the anomalous dimension of a given dia- 
gram plus the counter-terms and the insertion of the effervescent operators, 
i)-iii), the "complete" contribution. The advantage of combining different 
terms (bare diagram, counter-terms and insertion of effervescent operators) 
diagram by diagram is that this allows several checks on the contribution of 
any given two-loop diagram to the anomalous dimension. Thus for example, 
in HV, two diagrams which go one into another via a Fierz rearrangement, 
diagrams Vio and V12 in fig. 10 or P2 and F2 in fig.ll for instance, give the 
same "complete" contribution to 7*^^-* for left-left operators[|l^ . 

Similarly, relation (^TJ) can be shown to be true for the "complete" con- 
tribution of any single two-loop diagram. It is then possible to use it as 
a further check of the calculation in two different regularization schemes. 
Eg. (pT|) , used diagram by diagram, is a relation which connects single poles 
and finite terms at one loop with single poles at two loops. When satisfied, it 
authomatically ensures that the relation holds both for and aeOg correc- 
tions. A more extensive discussion of the checks done on the calculation of 
the single diagrams will be given in sections ^]4| and |5.5| . Before this we recall 
some basic facts about the HV and NDR regularization schemes, sec. [4.2| , and 
introduce the "effervescent" operators by considering the one- loop diagrams 
in these two regularizations, secs.|5.1[-|0|. Then we explain the construction 



of "complete" diagrams by two specific examples, one taken from the set 
of current-current diagrams, fig. 10, and the other from the set of penguin 
diagrams, fig.ll, sec.|5^ and |53. 
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4.2 HV and NDR Regularization Schemes 



In this subsection we report for completeness the computational rules of the 
HV and NDR regularization schemes, which have been used in the present 
work. In both the schemes Feynman diagrams are regularized by performing 
the integrals over the loop momenta in D = 4 — 2e dimensions. The two 
schemes differ in the definition of 75, which, in the case of the NDR regu- 
larization, can lead to inconsistencies and has to be treated with particular 
care. 

1. Naive Dimensional Regularization (NDR) 

In NDR 7-matrices are Z)-dimensional and only the D-dimensional met- 
ric tensor is introduced, following the convention: 

g,u = 9u, , <^?; = < , < = ^ , Tr(I)=4 (59) 

The D-dimensional Dirac matrices obey the usual algebra: 

{7^,74 = 2(7^, (60) 

and 75 anticommutes with all the 7's: 

{75,7m} = (61) 

The above definition of 75 may give problems in the evaluation of closed 
odd parity fermion loops, as for example Tr (75 7^ 7;^7p7o-), which are not 
unambiguously defined. For this reason it is not garanteed that it is 
possible to obtain the correct results by using the NDR scheme. In 
ref.0 and in the present calculation it is shown that, up to two loops, 
by fixing the ambiguity in the closed odd parity fermion loops, it is 
possible to obtain results which are in agreement with those obtained 
in the HV regularization scheme. In the present work we have fixed 
the ambiguity following the prescription of ref. We have verified 
that the prescription of ref.[Q gives the correct result for the triangle 



anomaly, as it is also the case for the HV scheme. 

2. t'Hooft-Veltman Regularization (HV) 

The HV regularization scheme has been proved to be unaffected by 
ambiguities or inconsistencies in the algebra of the 7 matrices ^7\, p^j . 
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The D-dimensional 7-matrices 7^^ are decomposed in two parts, a 4- 
dimensional part '^^ and a (D — 4)-dimensional part 7^^, i.e.: 

= + Y (62) 

with different anticommutation relations with respect to 75. The 4- 
dimensional 7's anticommute with 75, while the 7's commute with 75: 

{%,75} = , [7;„75] = (63) 



Besides the D-dimensional tensor of eg . (|59|) , two other metric tensors 
are introduced, namely the 4-dimensional tensor g^j^ and the (— 2e)- 
dimensional one g^y\ 

{lti,lu} = '^9tiu , {7m,%} = 2^m^ , {7M,7i/} = (64) 
with the following "mixed" contraction properties: 

r, = ^, 9', = -2e, c,% = Q (65) 

and 

« = ^;: , 9''x = K (66) 

Using HV for computing Feynman diagrams, it is convenient to take 
the external momenta in four dimensions (the loop momenta being 
integrated in D-dimension) . 

The bare operators inserted in the one- and two-loop diagrams are 
defined by using only 4-dimensional gamma matrices, i.e.: 

7^(1-75) ®%(l-75) or 7^^(1-75) ®%(l+75) (67) 

In the algebraic programs, written for the calculation of the diagrams, 
it is convenient to write the weak four-fermion operators in (|67D using 
(1 ± 75) as projectors: 

^(l+75)7'^(l-75) ® (1+75)7^(1-75) (68) 
and analogously for the 7^ ® 7^/? case. 
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In 4-dimensions, in the massless theory, all operators appearing in the 
weak effective Hamiltonian are suitable ffavour combinations of four-fermion 
left-left and left-right operators. This happens because long strings of 7 
matrices, appearing at one- or two- loops, can be reduced to left-left or left- 
right four-fermion operators, using the completeness of the Dirac algebra. 
This completeness is lost however using a dimensional rcgularization. In this 
case we can define an infinite number of independent four-fermion operators 
which, by quantum numbers, can mix with the original basis, at different 
orders in perturbation theory. The extra operators are those called previ- 
ously "effervescent" operators. They do not exist in four dimensions and are 
artefact of the rcgularization. It would be not correct however to compute 
the renormalization of the operators without taking into account the mixing 
of EO with the 4-dimensional operators. The reduction to the four dimen- 
sional basis can only be done at the end of the calculation of the two-loop 
anomalous dimension, as explained in the next section. 
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5 Diagrams and Counter-terms 



In this section we report the calculation of the one-loop diagrams in HV and 
NDR. This allows us to introduce the EO induced by the regularization and 
give the elements necessary to compute the coefficient functions C{Mw) as 
explained in sec.§. We also give some specific examples of the calculation 
of two-loop "complete" diagrams, both in the current-current and in the 
penguin case, in order to explain the method used in this work and the 
relations between the HV and NDR regularizations. A complete list of the 
contribution of all the diagrams is reported in Appendix. 



5.1 Current-Current Diagrams at One Loop 

Let us consider the diagrams in figs.3a-c, where the wavy lines can be due to 
a gluon or a photon exchange. The Dirac structure inserted in the vertex can 
be of the form 7^ ® 7^^ or 7^ ® 7^/5 and we consider both cases separately. 
After the loop integration, but before the simplification of the Dirac algebra, 
the pole terms of the diagrams have the following structure (for the two- 
loop anomalous dimension we only need to consider the pole contribution of 
4-dimensional and effervescent operators): 

G'-lLr = -C^^^ ^ {YYn ® lulpT,L,R + ® T,L,Rlul,) 

G^Ilr = ^^^'^^(7Vr^® W7p7. + r^7V®7.7pW) (69) 

where C^^\ C^"^^ and C^^^ are coefficients which summarize the colour/charge 
dependence of each diagram. and denote the weak vertex structure, 
given in NDR and in HV respectively by: 

^ _ / 7^(1-75), NDR , _j 7^(1+75), NDR 

'^-j 7^^(1-75), HV ^^-17^(1+^5), HV ^'^^ 



By reducing the 7-algebra in (p5|), we can separate the contributions of 4- 
dimensional operators from those of EO. By definition, the contributions of 
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EO correspond to those terms which vanish under a suitable projection on the 
four dimensional basis. There are several possible choices of the projection 
operators, which in general will define different renormalized operators, i.e. 
they will give different two-loop anomalous dimension matrices. The results 
can be easily related, by computing at one-loop the renormalized operators 
obtained by different projections and we have checked the consistency of the 
contribution of several two-loop diagrams obtained with different projections. 
For the sake of comparison, in the following we shall use the same projection 
as in refs.[||, |T9| : 



• for 7^ (g) — > Fll = 1r® IvR 

• for 7^ ® 7^K — > = (1-75) ® (1+75) 

Then, in order to project on 7^ ® liiL^ we take the following trace: 

Pll (7^ ® 7ml) = [7/^(l-75)7-(l+75)7^(l-75)7,(l+75))] (71) 

while, in order to project on 7^ ® 7^7?, we use: 

Pii? {ll ® 7mk) = [7A'(l_^5)(l_^5)^^(l+^5)(l+^5)] (72) 

In eq.(|7l|) the sum over v is intended in 4-dimension in HV and in D- 
dimensions in NDR. When projecting a string of 7-matrices, the values of the 
traces in (^) and (|72D are taken as normalization factors. With the projec- 
tors introduced in eq.(0) and eg. (|72D , we obtain the following decomposition 
for the one- loop vertex diagrams in (|69|): 



G^lIlr = C(^^^(i^SL/?(^)7i:®7.L,.? + 4lLii) (73) 

G^^Un = C(=^^^(i^SU(e)7i:®7.L,i^ + 4tL 
where the e-dependent coefficients are given, up to 0(e), by: 
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and: 

.(1) 



Fl'jiie) = 4^^^ h7 (76) 

p(2). . _ / -4(1 + e) NDR 

p(3). . _ / 16(1 -e) NDR 



are the contributions due to the EO. 



From eg. ([7^) , we obtain the counter-terms to left-left and left-right op- 
erators, including those proportional to EO: 



Gfl = C^^'^^(47^:®7.L + 4'2) (79) 

and similarly for the LR case. By subtracting the counter-terms in eq.(|79|) 
from the result of the calculation of the diagrams in figs.3a-c, we thus obtain 
the renormalized operator matrix elements expressed in terms of the tree- 
level matrix elements and of the matrix f (s), see eq.(^). The pole and 
finite terms coming from the calculation of diagrams Vi^a are reported in 
table H both in HV and NDR, for the insertion of a 7^ ® 7^l and 7^ ® 7^r 
operator. Only the terms proportional to the original operators are shown. 

Notice that the insertion of an operator 7^ does not produce "effer- 
vescent" operators when we compute the vertex renormalization of the vector 
current 7^, diagram Vi, since 75 is not involved. As a consequence, in HV, 
using the MS subtraction procedure, the finite term of the operator 7^ ® 7^ 
is different from the finite term of the operator 1/2^7^ ® 7^^ + 7^ ® 7/^^) • In 
particular for the vertex renormalization of 7^ one obtains: 

^(i + i)x7i;®7. (80) 

instead of: 

■ (- + ^) X ^ U ® + 7^ ® 1,r) (81) 



An^e 2' 2 
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diagram 


M 


7l ® 7a.l 


7l ® 7Mfl 






(lA) 


0{1)hv 


0{1)ndr 




0{1)hv 


0{1)ndr 


VI 


2 


1 


5/2 


1/2 


1 


5/2 


1/2 


V2 


2 


-4 


-9 


-9 


-1 


-3/2 


-7/2 




2 


1 


5/2 


1/2 


4 


10 


6 


Pi 


1 


-4/3 


-20/9 


-8/9 










1 


-4/3 


-20/9 


-20/9 


-4/3 


-20/9 


-20/9 



Table 2: Singular and finite terms for diagrams in figs.3a-c and 9, with a 
7l ® 7^1^ or a 7^ 7;^/? Dirac structure. The multiplicity of the diagrams 
is also reported in the table. Colour-charge factors and a the factor as/^'K 
(ae/47r) are omitted. 



as can be read in table ^ This difference has important consequences as will 
be discussed in the following. 

5.2 One-Gluon/Photon Penguin Diagrams at One Loop. 

The one-loop penguin diagrams are much simpler, due to the absence of EO 
at the one- loop level. The computation of the diagrams in fig. 9 gives, both 
in the HV and in the NDR schemes, the following Dirac structure (we omit 
the colour-charge factors): 

(Penguin) = - - ^ ^ {q'^ - Q'^l) (82) 

where is the momentum of the gluon/photon (see fig. 9) and is 4- 
dimensional or D-dimensional in the HV or in the NDR case respectively. 
We find: 

F^ie) = 1 + le 
Fr^'ie) = 1 + le (83) 
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We conclude that both in HV and NDR the counter-term has the following 
form: 

Gl, = ^\{q'n-q'i) (84) 

5.3 Two-Gluon or One-Gluon+ One-Photon Penguin 
One Loop Diagrams 

The presence of two gluon contributions is a typical feature of the two-loop 
calculation. The two-gluon (one gluon-|-one photon) diagrams, fig. 12c, only 
enter as counter-terms at the two-loop level, see for example fig. 15. These 
counter-terms exist only if, in any two-loop penguin-like diagram, we subtract 
completely the divergent part of the internal sub-diagram, without making 
use of the equations of motion. This is also true for the longitudinal compo- 
nent ~ q^(/i of the one gluon (photon) counter-term, which vanishes by the 
equations of motion. However one may choose to subtract only counter-terms 
which do not vanish by the equations of motion. The contribution of single 
diagrams will be modified, but the final result will be the same. Thus, for 
example, consider the diagrams P2 and P3. They can be computed with or 
without the two gluon counter-terms, Gf^ and G\g respectively, whose sum is 
equal to zero, fig. 16. The abelian part, i.e. the sum of the two "complete" di- 
agrams, remains the sameQ as schematically indicated in fig. 16, even though 
each of them is modified. Something similar happens to the longitudinal 
contribution of the penguin counter-term, eq. (P^ , for diagrams P14 and P15. 
If we subtract all the terms ~ 1/e in the penguin internal sub-diagram, we 
are also subtracting the longitudinal term q^^. We may however make use 
of the equations of motion and subtract only the term ~ Q^Il- The sum 
of the two diagrams (corresponding to the abelian case) remains the same 
because of the cancellation of the contributions of the counter-terms due to 
the longitudinal components, fig. 17. The same happens in the non-abelian 
case, but in a more complicated way which involves the difference of the two 
diagrams and also non-abelian diagrams. We have explicitely checked that, 
if we subtract only those counter-terms which do not vanish by the equations 
of motion, we obtain the same anomalous dimension matrix. However we will 
give the results by subtracting all the pole parts of the internal sub-diagrams 

That is their sum is the same with or without the subtraction of Gof". 
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because in this way eq. ( |5TD remains vahd diagram by diagram, see below. 

The cancellation of the sum of the counter-terms proportional to the lon- 
gitudinal contribution in eq. ( |MD and of the two-gluon (one gluon-one photon) 
counter-terms corresponds to the substitution p9|: 

7Ai7^^"i^^,.^2 ^ (^i7^^2) E ^7^^ (85) 

i.e. to the subtraction done using only four-fermion operators. Notice how- 
ever that minimal subtraction of the pole term proportional to the operator 
ipi'^lD^ F ^^jjil)2 in HV, is non-minimal in the basis where we use the 7^ ® "^^l 
and 7^®7^_R operators. This happens because, as explained above, the finite 
term of the renormalized operator 7^ ® 7^ is different from the finite term of 
the operator | {^1 ® 7^^ + 7^ ® 7^_r) . It is this last operator that we have 
indeed to subract for consistency with the current-current diagrams. The 
difference between 7^ ® 7^ and | [^1 ® 7^l + 7^ (S> 7^r) is important only 
for the counter-terms of diagrams Pig and Fig. Since we have to apply a 
non-minimal subtraction to the lower vertex, the contribution of Pig and Pig 
is non-zero, but it is given by one half of the pole of the one-loop penguin 
diagram times the difference of the finite terms of the two operators 7^ ® 7^ 

and \ (7l ® 7^^ + 7l ® 7^^) : 

1-4 41 

- 3^ X (-2, ^ -- (86) 

in agreement with eq. (^Uj) . This point was overlooked in ref.0. 

The results in Appendix have been obtained by subtracting the pole term 
of the internal divergent sub-diagrams, without making use of the equations 
of motion. With this choice, the relation between the anomalous dimension 
in NDR and HV remains valid diagram by diagram, eq. (|5l|) . We have checked 
that our results satisfy eq. ( pT] ) diagram by diagram, as the reader can verify 
by himself using the results from tables and |6|-|10|. Below we will discuss 
how this works in three specific examples. 

The fact that the sum of all the counter-terms due to operators which 
vanish by the equations of motion cancel, makes the relation (^) valid also 
for the complete matrix. 

For completeness we report also the explicit form of the two-gluon counter- 
term. The pole part, obtained from the first diagram in figs. 12c, is given by 
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(colour-charge factors omitted): 

4 F/^'^^^(e) 



g'^" (87) 



3 e 

Q,u ^ ^2)^ _ (2fc, + k^r^l + (fci + 2k2Y^l) (88) 

and the finite part, i.e. Fp^'^^^{e) is (as expected) equal to the finite part 
obtained from the diagram of fig. 12b, cf. eg. (p3|) . We thus find that the 
counter-term is in this case: 



G% = -Gi = ^Iq'^ (89) 



A diagramatic representation of the one-loop diagrams and the corresponding 
counter-terms is given in figs.l2a-c. 



5.4 Two-Loop Current-Current Diagrams. 

After the detailed study of the one loop diagrams, we are ready to show 
the construction of a "complete" two-loop diagram, including all necessary 
counter-terms. We will consider as an example the diagram of fig. 10 denoted 
as Vi7, where we insert a 7^ 7/^l vertex. We denote by -D^^ the double 
and single pole contribution from this diagram. We then substitute to the 
internal loop, including the string of Dirac matrices, the suitable counter- 
term, G^ll- We denote by C^^'^l double and single pole contribution 
from this counter-term. We have finally to correct for the mixing between 
4-dimensional operators and EO which occurs at one loop. This corresponds 
to the last term in eg. (^91) , the formula which gives the two-loop anomalous 
dimension[^ To obtain this term we substitute to the internal loop the 
combination : 

^(17) _ 1^(1) /^7Vr^p7.®r^LFS(o) 




4^ Fi'l 



(90) 



The last term in eq. (^) receives a contribution also from non " effervescent" operators 
in the case of a non minimal subtraction. In this respect there is not much difference 
between EO and finite subtractions of 4-dimensional operators. 
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The above equality is valid up to terms which do not contribute to the two- 



loop anomalous dimension because are of higher order in e, cf. eq.(79). 

(17) 

We denote by E(j £^ the contribution from this insertion. The contribution 
proportional to the first coefficient of the beta function (3q in eq. (^) is absent 
for the diagram V17. It will only give a contribution for those diagrams, 
which contain an internal loop corresponding to the renormalization of a^, 
as for example V12. The term ~ (3q is authomatically taken into account 
by subtracting the counter-term corresponding to the renormalization of the 
strong vertex. 

The result for the "complete" diagram is thus given by the sum of the 
three different contributions, see fig. 13: 

6i'? = z)£'-c<;,iU4"A (91) 

This ends our discussion of current-current counter-terms. In tables ^ and 
we report in units of (Q;s/47r)^, (astte/lSvr^), the double and single pole 
contribution of all the current-current diagrams and of the corresponding 
counter-terms, specifying the contribution of EO counter-terms whenever 
present. 

From table |^ we notice that, in HV, diagrams which can be changed 
one into the other by Fierz rearrangement, like Vn and V20 for example, 
give exactly the same contribution to the anomalous dimension (double and 
single poles). This is however true only for the "complete" diagrams, D\^i^ 
and D^^i\ while there is no relation between the corresponding bare diagram 
contributions, -D^x'* ^^cl D^^^\ This happens because the HV regulariza- 
tion scheme preserves the Fierz properties of the four dimensional basis of 
the renormalized operators defined via MS subtraction^ The renormahzed 
operators are however obtained only in combination with all the possible 
counter-terms. On the other hand, MS NDR does not respect the Fierz 
properties of the renormalized operators, as it can be seen from the results 
reported in tables It would be possible however to obtain a two-loop 

anomalous dimension which respects the Fierz symmetry by a suitable fi- 
nite one-loop redefinition of the renormalized operators [O, |l9i. Since Fierz 



^''This is strictly true with the projection operators (^Tj) and (JT^). Also in HV one can 
choose projection operators such that Fierz rearrangement relations are not satisfied at 
two loops. 
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rearrangement is not valid in NDR, 7^ ® 7^ij and —2(1 — 75) ® (1 + 75) 
Fierz rearranged operators renormalize differently. Notice that in the case 
of —2(1 — 75) ® (1 + 75) operators, diagrams denoted by Pi in fig.ll do not 
vanish. We have also computed the renormalization of the effective Hamil- 
tonian in the Fierz rearranged basis and verified that we obtain consistent 
results, i.e. that the anomalous dimension matrices in the basis (01) and in 



the rearranged one satisfy eq.(|5lD. 

As already anticipated, the relation which connects the two-loop anoma- 
lous dimension in two different regularization schemes (^TJ) holds diagram by 
diagram. This means that in eq.(pT|) (or in eq.(^)), we have to interpret 
the two-loop anomalous dimension 7(1) (or Zf)) as the contribution of the 
particular diagram under consideration. We have also to interpret the one- 
loop anomalous dimension 7*^°-' (or z[^^) and one- loop matrix f (or Zq^^) as 
due to all the possible diagrams obtained by eliminating one of the gluon 
propagators in the two-loop diagram. 

We illustrate further this point by discussing again the current-current 
diagram Vn of fig. 10. Eq. (|5ll) , for a single diagram, has to be interpreted as 
follows. The commutator Af,'jf'^ corresponds to the difference of the finite 
parts of the internal diagram Ary^, corresponding to diagram Vi, times the 
pole term of the external diagram minus the difference of the finite parts 
of the external diagram Ary, times pvi- Omitting colour-charge factors, in 
the case at hand we have: 

l-(-l) = ^((-2)x(-4)-lx0) (92) 
which coincides with the difference of the 1/e term reported in table ^ for 

In general the appropriate formula is: 

Ap(^) = P%n - V% = \ E [Ara ■ pT - " Ar,] (93) 

a, 6 

a and h indicate all the possible sub-diagrams which appear in the two-loop 
diagram. Pn^drhv^ which is the single pole of the two-loop diagram, has a 
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subscript because it depends on the regularization. p^^^ is the one-loop pole 
term for the diagram a. Ar^ is the difference of the finite terms between 
NDR and HV for the sub-diagram a. If p^^-* includes the poles corresponding 
to the renormalization of the strong coupling constant, eq.(^) incorporates 
authomatically the last term of eq. (|5ll) . We have verified diagram by diagram 
that all our results satisfy eq.(^). 



5.5 Two-Loop Penguin Diagrams. 

We will consider, as an example of the two-loop penguin diagrams, the dia- 
gram of fig. 11 denoted as P3, where we insert a 7^ '~i^L vertex. We denote 
by P2l double and single pole contribution from this diagram. We then 
substitute to the internal loop, including the string of Dirac matrices, the 
suitable counter-term, G)^^. In this case however we have also to subtract 
the two-gluon counter-term since the first of the two diagrams in fig. 12c is 
also contained in P3 as a divergent sub-diagram. We denote by Cp\j^ the 
double and single pole contribution from all the counter-terms. We have 
finally to correct for the mixing between 4-dimensional operators and EO 
which occurs at one loop. This corresponds to the last term in eq.(^), the 
formula which gives the two-loop anomalous dimension. As was the case for 
1^17, to obtain this term, we substitute to the internal loop the combination: 



£.(3) _ /r^iV^rvT^^SW (-4ra 



- -C^^^^xIeS (94) 

and denote by the double and single pole contribution from this inser- 

tion. We thus obtain: 

Pi'l = Pn - 41l + 41l (95) 

A diagramatic representation of the subtraction procedure is reported in 
fig. 14. In HV two diagrams which go one into the other by Fierz rearrange- 
ment give the same contribution to the anomalous dimension. All the results 
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for the "complete" penguin diagrams are given in tables p|-p!0| for HV and for 
NDR. 

We can check eq.(^) also in the case of P3. We call Ar2g and p2g the 
finite part and pole term of the two-gluon diagram of fig. 12c and ^r2gct and 
P2gct the corresponding quantities for the insertion of the operator Q^^ of 
eq. (|88|) , see fig. 15. One then obtains {j)2g = —4/3, p2gct = 1, ^^^29 = 4/3 and 
^r2gct = 0): 

Ap(i) = ^ (Arvjpl - Arp, + Ar2gP2gct - P2gAr2gct) = y (96) 

in agreement with the results of table 

We notice that our results for the diagrams Pg and Fg in HV satisfy 
eg . (|5T|) , as can be easily verified. For these diagrams, in HV, the authors of 
ref.0 found a different result which fails to satisfy eq. (|51|) . 



5.6 Terms Which Vanish by Equations of Motion and 
Gauge Invariance 

The results reported in tables |6|-p!0|have been obtained by subtracting the pole 
term, including the EO, eqs.([79D, for all the divergent sub-diagrams. In doing 
so, we have subtracted also operators which vanish by the equations of motion 
or non-gauge invariant operators. As already discussed before, we could have 
subtracted only four-fermion operators and obtained the same final result. 
However contributions of single diagrams would have been different and the 
check between HV and NDR, eq. (pT]) , would have not be valid diagram by 
diagram. Besides this we have of course to take into account the mixing with 
"effervescent" operators. This mixing is responsible, besides other effects, of 
the non-vaninshing of diagrams Pie and Pie in HV. 

The presence of operators which vanish by the equations of motion and 
non-gauge invariant operators is signaled by the appearence of various tensor 
products in the calculation of the different two-loop Feynman diagrams (we 
follow the notation of ref.0): 

Ti = 7^®7'^ 
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When we take into account the multiphcity. 
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T2 = =^(1-75)®^^ 

Ts = Tix^ (97) 

T4 = {H'^-l5)^f + f{l-l5)®^)\ 



where —p is the ingoing momentum, r the outgoing momentum, q = p + r 

and Sf,au = Itilalu - l^lalti- 

T2 and T4 vanish by using the equations of motion: 

u{r)f = ^u{p) = (98) 

and 

T3 y T5 ^ -Ti (99) 
Non-gauge invariant operators can be ehminated by using the background 



field Feynman gauge ||3(| for the non-abehan penguin diagrams, for example 
P4 and Pq in fig. 11. This allows simple checks of the gauge invariance of the 
final result, because non-gauge invariant operators cancel when combining 
together sub-sets of two-loop penguin diagrams. 

We give an explicit example of this cancellation by considering diagrams 
P2, P3 and P4 in HV. The results for these diagrams, using the background 
field gauge, putting = = and reporting only double and single poles 
are: 

-|- (p r) 



36 




+ {p ^ r) 




+ {p ^ r) 



The explicit results are given for a gluon propagator stemming from the 
upper incoming quark leg, as shown in fig. 11. p is replaced by r (and the 
term proportional to e^^pa^ changes sign) when the gluon is attached to the 
outgoing quark leg. Notice that g^-P2' is different from zero and similarly 
for and P^, indicating the presence of non-gauge invariant operators^ 
However: 

q, {P^ + P-^) = 
q,{P^-P^ + 2P^)=Q 

if we use the backgroud field gauge for the non-abelian diagram. The sum 
of the above diagrams eliminates the non-gauge invariant operators. By 
writing the term ~ e^^pa as a term proportional to g^7j^ — (/iLq^ plus terms 
which vanish by the equations of motion and by eliminating for the same 
reason all terms proportional to |^ or we arrive to the results reported in 
the tables. 

We conclude this section by summarizing some checks done to verify the 
correctness of our calculation. 

1) We have verified that the 1/e^ contribution of the counter-term is twice 
the corresponding term of the original diagram, as imposed by eq.(^7|). 

""^^We use the on-shell identity p ■ q — /2 
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2) We have verified the cancellation of all the single poles accompanied by 
logarithms of the external momenta, i.e. ~ 1/e x /n(p^//i^), 1/e x /n(r^//i^) 
or 1/e X ln{q^/ ^^). For penguin diagrams, this cancellation does not follow 
authomatically from 1). 

3) We have verified the cancellation of all non-gauge invariant terms in 
the background Feynman gauge and repeated the calculation in the standard 
Feynman gauge with the same final result. 

4) We have verified that we get the same result by subtracting all the pole 
terms in the internal loops, i.e. by subtracting also operators which vanish 
by the equations of motion, or by subtracting only those counter-terms which 
are proportional to four-fermion operators. 

5) We have verified that eg. , which relates the NDR and HV schemes, 
is valid diagram by diagram. We have also verified that eq.(51) is valid for 
the matrices 7^^^ and 7^^^. 

In doing so we have found a difference in the result of diagram Pg (Fg) 
between our calculation and the calculation reported in ref.0], as can be 
read from tables P-|TO|. Our results for these diagrams, Pg and Pg, satisfy 
eq. (|5T|) , while the results of ref. p| do not. Moreover the authors of ref.p 
overlooked the different renormalization of the operator 7^ ® 7^ with respect 
to the renormalization of 1/2(^1 ® 7^/, + 7^ ® 7^_r) and consequently do not 
consider the diagrams that we call Pig and Pig. These differences compensate 
in the calculation of 7^-*. It is however easy to show that, by using for the 
diagrams the results of ref.[@, one cannot satisfy the relation (|5l|) for 7^^''. 
As a further check we also computed the anomalous dimension matrix at 
0{a1), since in this case different diagrams enter with different weights with 
respect to 7^^^ and 7^"^^. We found again that only our results satisfy (pT|), 
contrary to those of ref.0. This gives us further confidence on the results 
reported here. 

We finally notice that in ref. the authors have only computed the 
anomalous dimension matrix in NDR. They get the matrix in HV using 
eq. (|5T|) . Since we agree in NDR, the final result of ref. is correct, in spite 
of the difference found for the two diagrams in HV. 
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5.7 Quark Self-Energy Diagrams 

The calculation of the quark self energy diagrams at one and two loops is 
straightforward and has already been done by several authors, see for example 
IS]. We only report in table ^ of the Appendix the results of the different 
diagrams shown in figs. 4, 6 and 18. 
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6 The anomalous dimension matrices 



In this section we collect the results of the calculation of the anomalous di- 
mension matrices at one- and two-loop level both in the HV and NDR renor- 
malization schemes and the matrices Af and As, relevant for the comparison 
of the results in the two schemes. 



6.1 One- loop results 



The one-loop anomalous dimension matrices |T5|, and Af are presented 
here. They can be computed from the pole and finite parts of the one-loop 
diagrams in figs.3a-c and 9, given in table ^ and from the diagrams of figs. 4 
and 6a, table ^. We also report the 6f and 6s matrices introduced at the end 
of sec.|] to take into account the bottom threshold. 

The one-loop 0{as) matrix is given by 



7^ = (100) 



6 


6 


Af 


6 


6 


N 




























































2 


2 


2 


22 


I 


3N 
4 


m 

18N~2f 
3N 


3 

-18+2Nf 
3N 


3N 








6 

N 


-— f 

3NJ 


-f 













2u~d 


2u-d 


2u-d 








lu—d 


3 

2u-d 


lu—d 


3N 
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3N 


















2 














3 














3'' 














-6 














18Af^+2Af/ 
3Af 






6 

N 




-6 










2u-d 

\ 
3 






N 






6 

6 



6 


2u-d 
3 








6 

Af 



v 

where /, u and d represent the number of flavours, the number of charge 2/3 
up-like quarks and the number of charge —1/3 down-like quarks respectively. 

For the electro-magnetic anomalous dimension matrix one finds: 
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7, 



(0) 



(101) 





8 
3 






















8 
3 












































































































4 

3 






















4 
3 










4 
3 











[ 











4 
3 









16Af 

I 
, 27 

-16+(16n-8d)A^ 
27 

-16N+16u~8d 

, 27 
72+(16M-8d)Af 

27 
16M-8d 

27 

36+(16M+4d)Af 

27 
16M+4d 

. 27 , 
8+(16M+4tj)Af 

27 

8Af+16M+4d 
27 









8 
3 



4 
3 







i 

27 

-88+(16M-8d)jV 
27 

-16jV+16M-8d 
27 

(16n~8d)Af 

27 
16M-8d 

27 

(16n+4d)Ar 

27 
16M+4d 

. 27 

-28+(16M+4d)Ar 
27 

8Af+16M+4d 
27 



At one loop, from the finite terms in table ^ we also obtain: 
Af = 



\ 




_8 

3 







-1/ 



V 



-2Ar2 
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-2 


























-2 



4-27V2 
N 
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3Af„ 
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3 „ 
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(102) 



in the strong case and 
As = 



(103) 



41 





_2 

Q 






















2 
9 






















2 

3 






















2 
3 






















2 
3 






















2 
3 
















10 
9 






















10 
9 










2 
9 











V 











2 
9 









87V 

27 
_8_ 

27 

8N 

27 
20 

9 



4 
9 



4_ 

ik 

27 









20 
9 



4 

9 






SN 

27 
_8_ 

I 
27 

SN 

27 








\ 



4 
9 









4f , 

27 9 / 



in the electromagnetic one. 



5f and 6s in 


eq.® 
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where u (d) is the number of active u-type (d-type) quarks. These matrices 
are equal in the HV and NDR schemes. 



6.2 Two-loop results in HV and NDR 

All the ingredients needed for the calculation of the two loop anomalous 
dimension have been given in sees. ^ and ^ Our final results for the two- 
loop anomalous dimension matrix introduced in eq. (0) are given below in the 
HV and NDR regularization schemes. 

Each Feynman diagram can be schematically represented as the product 
of three factors: 

(Colour — Charge) x (Dirac) x (Numerical result of loop integrations) 

The numerical contribution of each diagram can be found in tables for 
vertex and penguin diagrams. It is given in units of as/47r (ae/47r). The 
Dirac structure can be read directly from the diagrams with few exceptions 
which are explicitly reported in the tables. The results of the easy but tedious 
calculation of the colour-charge factors are not given. 

The non-zero two-loop anomalous dimension of the weak current has been 
put to zero by a suitable redefinition of the current at one-loop. The effect of 
this redefinition has been already taken into account in in the calculation of 
the coefficient functions. Without this redefinition, using the minimal MS 
subtraction in the HV scheme, we would have found 7j = (A^^ — 1)/2N x 
4/3o. Both the one- loop and the two- loop anomalous dimension matrices are 
comprehensive of the self-energy contribution. 
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Table 3: Elements of the two- loop anomalous dimension matrix (7^^'),^ . 
elements which are not reported are equal to zero. (Continue) 
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Table 3: (Continued) Elements of the two- loop anomalous dimension matrix 

('Ts ) ij ' 
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1,1 

1,2 
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25 

3 
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2,2 
2,3 
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2,5 
2,6 
2,7 
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2,9 
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556 



3Ar 



243Af 
268 
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268 



243JV 
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520Ar _|_ 884 

243 243Af 
_ 52 
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520 Af _ 1636 
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556 
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27 
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(u — 
u — 
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56 



+ 
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243 



243Ar 243 V 
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16 
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1688A f 
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1690 
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27 + ~ ~ 2^ 
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32 
27 
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Table 4: Elements of the two- loop anomalous dimension matrix (7^ ^)jj- The 
elements which are not shown are equal to zero. (Continue) 
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32 



27-' ' 243 

J 243N " 
340 „ 



100 



243Af 
100 



243 
^ d 
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Table 4: (Continued) Elements of the two-loop anomalous dimension matrix 
(7^^^).. (Continue) 
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\d 
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\d 
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243 

U - 
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28 _|_ 88Af 
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\d 
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\d 
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232A' 
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1100 AT 
243 



794 
243 AT 
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_32 I 88iV 
27 27 
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-§-f 



740Af 
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38 
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Table 4: (Continued) Elements of the two-loop anomalous dimension matrix 
(7f )),,. (Continue) 
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ihJ) 


HV 


NDR 


(10,3) 

(10,4) 
(10,5) 
(10,6) 
(10,7) 

(10,8) 
(10,9) 

(10,10) 


243 2437V ' 243Af 
287iV 20 1 268 14 j 

243 3N 243 243 
44 268 ^ I 14 ^ 
243 243iV ' 243JV 
44Ar , 268 14 j 

243 ~ 243 243 
AN^ 8 _|_ 520W 1 

Q 27 243 

884: 64iV J 1 260 j 

2A3N"' ' 243 ' 243Af 
_4iV _ ^ 4 . 

27 9 " 3" 

4Af2 _|_ 8 _|_ 520Af 

9 3 243 
1636 1 64Af J 316 j 

243N"' ' 243 243Af 
22Ar 2 I 124 I 28 J 

9 3Ar ~ 27 " ~ 27" 


1000 000 ^ ID ^ 

243 243Af 243Af " 
287iV 20 I 388 1 16 ^ 

243 3N 243" 243" 
44 388 16 J 

243 243JV 2437V " 
44Ar , 388 1 16 J 

243 ~ 243 ^ 243 
AN^ 32 1 280iV _|_ 

9 27 243 ~'~ 
620 32Nj 40 J 

243Af 243 243Ar 
20Af lOD,, I 8 ^ 

27 27 ~ 27 

4Af2 _|_ 32 I 280Af 

9 ~'~ 27 243 
1612 32Ar J 328 j 

2A3N " 243 243Af " 
26Ar 2 I 148„, I 40 J 

27 3Ar ~ 27 " ~ 27" 



Table 4: (Continued) Elements of the two-loop anomalous dimension matrix 
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Figure Captions 

• fig.l: Tree- level current-current diagram. 

• fig.2a-c: One-loop 0{as) current-current diagrams in the full theory. 

• fig.3a-c: One-loop 0{as) current- current diagrams in the effective the- 
ory. In the figure the relevant momenta of the external legs are given. 

• fig.4: One-loop 0{as) quark self-energy diagram. 

• fig.5a-d: One loop 0{ae) current- current diagrams in the full theory, in- 
cluding the non-abelian diagrams. The Z^-W diagrams are also shown. 

• fig.6a-b: The quark and W self-energy diagrams at 0{ae). 

• fig. 7: Box diagrams contributing to B{xt). 

• fig.8a-b: Penguin diagrams in the full theory. In the electro- weak case 
also the non-abelian diagram is shown. 

• fig.9: Penguin diagrams in the effective theory. In the figure the rele- 
vant momenta of the external legs are given. 

• fig. 10: Current-current diagrams at two loops. 

• fig. 11: Penguin diagrams at two loops. 

• fig. 12: Schematic representation of the counter- terms: current-current 
(a), one-gluon penguin (b) and two-gluon counterterms (c). 
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fig. 13: Diagramatic representation of tlie subtraction procedure for a 
current-current diagram. The "complete" diagram is obtained by sum- 
ming the diagram with a bare operator inserted, the counter-terms 
including those corresponding to effervescent operators and the term 
defined as E^j^p in eg . (pOD . 

fig.l4:Diagramatic representation of the subtraction procedure for a 
penguin diagram. The "complete" diagram is obtained by summing 
the diagram with a bare operator inserted, the counter-terms including 
those corresponding to effervescent operators and the term defined as 
E^l in eg. (P^ . Also the two-gluon counter-term is shown. 

fig. 15: One loop diagram for the s ^ d + g + g (g + 'j) operator. From 
this diagram one can compute p2gct and Ar2gct, cf. eg. (PBD . 

fig. 16: This figure shows that the sum of diagrams P2 and P3 is the 
same with or without the two-gluon counter-terms. This corresponds 
to the cancellation of counter-terms which vanish by the eguations of 
motion in the abelian case. 

fig. 17: The contribution of the longitudinal term ~ q^^ of the penguin 
counter-term cancels when we sum the diagrams shown in this figure. 
This corresponds to the cancellation of counter-terms which vanish by 
the eguations of motion in the abelian case. 

fig. 18: Quark self-energy diagrams at two loops. 
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Appendix 



Diagram 


1 
e 


0(1) 


5*0 


-1 


-1/2 


^1 


1/8 




S2 


//2 




S3 


5/8 






1/4 






11/4 





Table 5: Single pole and finite part for the one-loop and two-loop self-energy 
diagrams in fig. 4, 6 and 18. For the two-loop case only the pole part is 
given. 
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Table 6: Two-loop pole contributions for the 7^ ® four-quark type dia- 
grams in fig. 10. For V29, V30 and V31, the results proportional to the number 
of colour N (V^) and flavour / (Vf) are separately reported. 
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Table 7: Two-loop pole contributions for the 7^ four-quark type dia- 
grams in fig. 10. For V29, V30 and V31, the results proportional to the number 
of colour N (V^) and flavour / (Vf) are separately reported. 
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Table 8: Two-loop pole contributions for the P-type penguin diagrams in fig. 
11 when a 7^ (g) 7^^ structure is inserted in the upper vertex. All penguin 
diagrams have a 7^(8)7^ structure, except Pio, Pu, Pu and P15, for which we 
explicitly give the 7^® 7^ (P^^) and the 7^(8)7^75 (P^"^) part. For Pi and Pr 
both the Feynman (P) and the background Feynman (P''^) gauge results are 
reported. Diagrams P5, P7 and P12 are not included, because their pole parts 
vanish and they do not contribute to the two-loop anomalous dimension. 
Some contributions to diagram P13 are identical for the bare diagram and 
the counter-term. This means that one does not need to compute these 
terms since they cancel in the final result. For this reason they have not 
been reported in the table. 
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Table 9: Two-loop pole contributions for the F-type penguin diagrams in fig. 
11 when a 7^ (g) 7^^ structure is inserted in the upper vertex. All penguin 
diagrams have a 7^(8)7^ structure, except Fiq, Fn, F14 and F15, for which we 
explicitly give the 7^®7^ (F^^) and the 7^(8)7^75 (F^^) part. For F.^ and F^ 
both the Feynman (F) and the background Feynman [F'^^) gauge results are 
reported. Diagrams F5, F7 and F12 are not included, because their pole parts 
vanish and they do not contribute to the two-loop anomalous dimension. 
Some contributions to diagram F13 are identical for the bare diagram and 
the counter-term. This means that one does not need to compute these 
terms since they cancel in the final result. For this reason they have not 
been reported in the table. 
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Table 10: Two-loop pole contributions for the F-type penguin diagrams in 
fig. 11 when a 7^(8)7/^^? structure is inserted in the upper vertex. All penguin 
diagrams have a 7^(8)7^ structure, except Fio, Fn, -F14 and F15, for which we 
explicitly give the 7i®7^ (F^^) and the 7^(8>7/i75 (-F'^^) part. For F\ and 
both the Feynman (F) and the background Feynman (F''^) gauge results are 
reported. Diagrams F5, F7 and F12 are not included, because their pole parts 
vanish and they do not contribute to the two-loop anomalous dimension. 
Some contributions to diagram F13 are identical for the bare diagram and 
the countcr-tcrm. This means that one does not need to compute these 
terms since they cancel in the final result. For this reason they have not 
been reported in the table. 
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